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Based on low temperature resistivity, heat capacity and magnetization investigations we show that 
the unusually strong suppression of superconductivity in Lu x Zri_ x Bi 2 BSC-type superconductors 
in the range x<0.08 is caused by the emergence of static spin polarization in the vicinity of non¬ 
magnetic lutetium impurities. The analysis of received results points to a formation of static 
magnetic moments with per Lu-ion. The size of these spin polarized nanodomains was 

estimated to be about 5 A. 

PACS numbers: 74.70.Ad, 74.62.-c 


1. The discovery of superconductivity at Tc«39 K in 
MgB 2 [1] stimulated a significant interest into the studies 
of a wide class of the alkaline-earth, rare-earth and 
transition-metal borides. Among them, in the family of 
higher borides RB 12 , zirconium dodecaboride (ZrB 12 ) is a 
BCS superconductor with the highest Tc ~6 K [2, 3]. An 
intriguing detail established for ZrBi 2 is the formation 
of Cooper pairs through quasi-local vibrations involving 
Zr 4+ ions located within truncated JJ 2 4 octahedrons in 
the UB 12 -type fee crystal structure [2]-[5]. In studies 
of the Einstein phonon mediated superconductivity in 
ZrBi 2 authors of [2]-[9] argue that s-wave pairing is 
characteristic for this compound, and that in this case the 
Ginzburg-Landau parameter k is located in the nearest 
vicinity of the threshold value nc= 2 -1 / 2 . Moreover, 
a crossover from type-I to type-II /1 superconductivity 
with temperature lowering was deduced in [3] from heat 
capacity and magnetization measurements. In contrast, 
in [10] the superconductivity in ZrBi 2 was interpreted in 
terms of d -wave pairing and a two-gap type-II regime 
was identified with parameters k p = 3.8 and Kd— 5.8. 
Additionally, a large size pseudo-gap (A~7.3 meV) has 
been detected employing high resolution photoemission 
spectroscopy in ZrBi 2 above Tc , and proximity to the 
quantum fluctuation regime was predicted from ab initio 
band structure calculations [11]. Thus, also similarities 
with cuprate high temperature superconductors (HTSC) 
may be supposed to revive the interest into studies of this 
low temperature superconductor. 

In case of non-magnetic impurity substitutions and 
their impact on superconducting properties one can 
choose between various scenarios, because for these 
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defects the pair-breaking mechanism in various models 
differs. According to Anderson’s theorem [ 12 , 13] 
and its extension to non-s-wave superconductivity (see 
e.g. [14, 15]) already a small amount of non-magnetic 
impurities can dramatically suppress superconductivity 
by pair-breaking in the case of anisotropic gap in 
a d -wave superconductor. Moreover, experiments on 
cuprates reveal [15] that a spinless impurity (Zn, Li etc.) 
introduced into a HTSC host produces in its vicinity 
a large and spatially extended alternating magnetic 
polarization. On the basis of NMR and fi SR spectra 
it has been demonstrated that this impurity-induced 
magnetization on the nearest neighbor Cu atoms in 
cuprates is associated with a dynamic moment [15]. 

In the case of RBi 2 the replacement of non-magnetic 
ions of Zr by Lu produces an about 15 times’ reduction 
of superconducting transition temperature (Tc~0A K 
for LuBi 2 [4, 9, 16]), and the origin of this large Tc 
suppression is not cleared up to now for these two 
compounds with similar conduction bands and crystalline 
structures. Indeed, inelastic neutron scattering studies of 
the phonon spectra in LuBi 2 and ZrBi 2 [5] have detected 
noticeable, but not dramatic changes in the position 
of almost dispersion-less quasi-local mode (15 meV and 
17.5 meV, correspondingly), which was proposed to 
be responsible for Cooper pairing. Only a moderate 
difference in the electronic density of states of these 
two compounds is caused by filling of a wide enough 
conduction band (~1.6-4-2 eV) when Lu 3+ — ion is 
changed to Zr 4+ in the RBi 2 unit cell, resulting to an 
elevation by about 0.3-0.4 eV of the Fermi level Ef for 
ZrBi 2 in comparison to LuBi 2 [11, 17]. 

Here we probed the evolution of superconducting 
transition temperature Tq and the normal state 
parameters for substitutional solid solutions 
Lu^Zn^Bia employing resistivity, heat capacity 
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Fig. 1: (Colour on-line) Temperature dependences of 
(a) resistivity, (b) magnetization at H =5 Oe and 
(c) specific heat in vicinity of Tq, and (d) the 
suppression of superconductivity Tq(x) in Lu^Zri-^B^. 


and magnetization measurements. It will be shown 
that the non-magnetic Lu impurity substitution (having 
a 4/ 14 configuration) produces a strong static spin 
polarization in the RB 12 matrix in the vicinity of 
lutetium ions. Simultaneously with the emergence of 
static magnetic moments with a value of about 3 /xb per 
Lu ion the received experimental results exhibit a strong 
suppression of superconductivity in Lu a ,Zri_ a; Bi 2 . 

2 . Studies of resistivity, transverse magnetoresistance, 
heat capacity and magnetization of high-quality single 
crystals of Lu a ,Zri_ a; Bi 2 solid solutions with x<0.08 
were carried out at temperatures in the range 1.8- 
300 K, in magnetic fields of up to 90 kOe (H||(001)). 
A standard dc four probe technique was applied for 
resistivity investigation with the orientation of measuring 
current I|| (110). Magnetization and heat capacity were 
measured using a PPMS-9 (Quantum Design). The 
single crystals of Lu a ,Zri_ a; Bi 2 were grown by vertical 
crucible-free inductive floating zone melting in an inert 
gas atmosphere. To verify the samples’ quality and 
the Lu content X-rays Laue back-pattern (see Fig.SI 
in the Supplementary Information) and microanalysis 
techniques were applied. 

3. Temperature dependences of resistivity (fig. la), 
magnetization (fig. lb) and heat capacity (fig.lc) show 
superconducting phase transitions with Tq in the range 
4.5-6 K for Lu 2 ,Zri_ a: Bi 2 (x<0.08) solid solutions. The 
resistivity p(T) drop to zero below Tc is accompanied 
both with the appearance of Meissner state diamagnetic 
response on M(T) curves and the stepwise changes in 
the specific heat C(T) (see fig.l and figs.2a-b, H =0 
curves). With the increase of lutetium content both 
the residual resistivity and the normal state specific heat 
rise dramatically (by a factor of 4, see fig. la and figs. 2a 
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Fig. 2: (Colour on-line) Temperature dependences 
of specific heat for Lu x Zri_a;Bi 2 with (a) x=0.007 
and (b) x=0.074 in magnetic field up to 90 kOe. 
Panels (c) and (d) demonstrate the separated magnetic 
component C m of heat capacity for Lu content 
x=0.007 and x=0.074, correspondingly. The fitting 
of Cm (T,H 0 ) by three-level Schottky relation (see 
scheme on panel (c)) is shown by solid lines. The 
splitting energies AEi(H), A E 2 (H) and concentrations 
A T TLSm(H) are shown for Lu x Zri_a;Bi2 with x=0.007 
and 0.074 on panels (e),(f) and (g),(h), correspondingly. 


and 2b), and the later one demonstrates a combination 
of superconducting step-like and Schottky-type anomaly 
(see e.g. fig. 2b, H =0 curve for x=0.074). It is worth to 
note that the Schottky anomaly for Luo.o 74 Zr 0 .o 926 E>i 2 
(fig. 2b, H= 0) is very similar to that one observed 
previously for LuBi 2 [9, 18, 19], and it may be interpreted 
in terms of formation of two-level systems (TLS) in the 
disordered RB 12 matrix of this cage-glass compound. 
The Tc(x) dependence for Lu a; Zri_ a ,Bi 2 solid solutions 
is summarized in fig. Id. It should be stressed that the 
suppression of superconductivity by Lu substitution is 
unusually strong (^0.21 K/at.% Lu), when supposing 
the doping of an s-wave superconductor by non-magnetic 
impurities. 

In external magnetic field the amplitude of the low 
temperature Schottky anomaly increases essentially and 
C(T) maximum moves to higher temperatures (for 
example, figs. 2a, 2b show the data for x=0.007 and 0.074, 
correspondingly). Evidently, the Zeeman component 
is a magnetic contribution to heat capacity and it 
may be separated from vibrational and electron heat 
capacity terms following the approach developed in 
[18, 19]. The resulting C m (T,H 0 ) dependences are 
presented in figs. 2c,2d for these two Lu contents together 
with their approximation (see solid lines in figs. 2c,2d) 
based on a three-level Schottky relation (see eq.(Sl) in 
supplementary materials) (the scheme presented in insert 
to fig. 2c). Both splitting energies A Ei(H), A E 2 (H) and 
the concentration of magnetic Schottky sites NTLSm(H) 
as deduced from the approximation are shown on figs.2e, 
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Fig. 3: (Colour on-line) Dependences of magnetization 
vs. temperature at H =2 kOe (a) and vs. magnetic 
field at T= 2 K (b) for Lu x Zri_ x Bi2. Approximations 
by Curie-Weiss relation (1) (panel (a)) and by Brillouin 
dependence (2) (panel (b)) are shown by solid lines. 


2g and figs.2f, 2h, correspondingly. The slope of straight 
lines on figs.2e and 2f (in the range H <40 kOe) allows 
to determine the g-factors from relation AEi^/isgiH/kB 
and obtain the values gi~3.2 and g2~7.8 for both studied 
crystals. Moreover, the concentration of TLSm is found 
to be within experimental accuracy equal to x values in 
Lu a ,Zri_ a; Bi2. From these results it may be concluded 
that magnetic sites are created by Lu substitution. 

The magnetic response in the normal state of 
Lu a ,Zr 1 _ a; B 1 2 superconductors was investigated by 
magnetization M(H,T) studies. Figs. 3a and 3b 

demonstrate the temperature dependences M(H 0 ,T) 
recorded at H o=2 kOe and the magnetization vs. 
magnetic field curves M(H,T 0 ) measured at Tq=2 K, 
correspondingly. It can be seen from fig. 3a, that a small 
(~4xlCC 4 ps/Zr) and about temperature independent 
Pauli-like paramagnetic response, which is typical for 
ZrBi2, changes into a Curie-Weiss-type magnetic signal 
originating from localized magnetic moments induced 
by the Lu substitution. Fitting the temperature 
dependences by Curie-Weiss relation 


M= Xo H + ^ eff H/(k B (T-e)) (1) 

where xo - temperature independent susceptibility, ks 
- Boltzman constant and 0 - Curie-Weiss temperature 
(see solid lines in fig. 3a) allows to determine the Curie 
constant and subsequently to estimate the effective 
magnetic moment p, e //~3.3-3.8pB per Lu ion (see 
also Tab.SI in the Supplementary materials). These 
values correlate very well with the saturated moments 
Ps~3-4/xb per Lu ion obtained from the analysis of 
magnetization vs. magnetic field dependences (fig. 3b). 

Moreover, in the regime of isolated magnetic impurity 
for £=0.007 the magnetization may be described with a 
good accuracy by the Brillouin-type dependence 


Fig. 4: (Colour on-line) Magnetic field dependences 
of normal state magnetoresistance Ap/p(H,T 0 ) 

(a) for Lu a; Zri_ x Bi2 with various Lu content at 
To=4.5 K, and (b) for £=0.015 at temperatures 
in the range 1.85-8 K. Approximations in the 
framework of (eq.4) are shown by solid lines in panel 

(b) ; curves are shifted by 0.2% for convenience. 
Panels (c) and (d) demonstrate temperature 
dependences of A and go parameters in (eq.4). 


M = N m g^ B J Bj (g/j, B JH/k B T) (2) 

(where N m —concentration of magnetic sites, 
Bj— Brillouin function, ps—Bohr magneton and 
J—angular momentum) with A r m ?»0.48£, gss3 and J~2. 
As a result, it encourages us to suggest as a simplest 
scenario that impurity ions of Lu are combined into 
dimmers with a magnetic moment ps=gJ=6ps and with 
a size equal to about d(Lu-Lu)~5.3 A — the distance 
between heavy ions in the RB12 matrix, which can be 
considered as a formation of spin-polarized nanodomains 
in Lua;Zri_ x Bi2. 

Magnetoresistance A p/p experiments were also 
performed to characterize the magnetic moments in 
these compounds with metallic conduction. According 
to Yosida calculations carried out within the framework 
of s-d exchange model, an appearance of negative 
magnetoresistance (nMR) is expected in the regime of 
charge carrier scattering on localized magnetic moments 
[20]. Both field and temperature dependences of nMR 
are controlled by local magnetization Mi oc through the 
relation 

- A p/p ~ Ml c . (3) 

The nMR effect may be considered as an independent 
argument in favor of local moments’ formation in 
a metallic matrix. Fig. 4 shows the results of 
magnetoresisance measurements on Lu x Zri_ x Bi2. In 
the normal state of ZrB 12 the main contribution to 
A p/p(H,To) is positive and it can be described with a 
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good accuracy by the well-known relation A p/p~g,D 1 2 3 4 5 6 7 8 9 10 11 12 H 2 , 
where p, B is the carriers’ drift mobility. Both the 
substitutional disorder and the cage-glass effect in 
Lua;Zri_ x Bi2 decrease dramatically both the mobility 
and amplitude of the positive component when Lu 
concentration increases (fig. 4a). For x>0.015 the 
appearance of nMR contribution to A p/p becomes 
evident, and for x~0.07 the negative term in applied 
magnetic fields (fig. 4a) prevails. Fig. 4b presents a set 
of magnetic field dependences A p/p(H,T 0 ) obtained for 
x=0.015 at temperatures between 1.8 K and 8 K. In 
the range between 8 K and 4 K a crossover from 
positive magnetoresistance to a combination of positive 
and negative components is observed, and similar to 
the approach developed in [21], these data may be 
approximated very well by relation 

A p/p = ~{Agp B JBj{gp B JH/k B T)) 2 + p 2 D H 2 . (4) 

Fitting data by relation (4) (see solid lines in fig. 4b) 
allows to estimate the g-factor g«3.5 and the angular 
momentum J=2. The obtained values are very close to 
those deduced above from the analysis of magnetization 
results by relation (2). Additionally, we have evaluated 
the behavior of coefficients A(T) and pd(T) in (4) 
(see figs. 4c and 4d, correspondingly). As can be seen 
from fig. 4c, the scattering of charge carriers on localized 
moments in vicinity of impurity sites appears below 8 
K and increases drastically with temperature lowering. 
On the contrary, the Pd(T) dependence demonstrates 
only a moderate elevation with temperature decrease 
(fig.4d). Thus, these positive and negative components 


of magnetoresistance become comparable at helium 
temperatures and in magnetic fields below 40 kOe. 

When discussing the possible scenario of the formation 
of spin-polarized nanodomains in this case, it is worth to 
note the mechanism of spin-polaron formation proposed 
for the rare earth higher borides with a cage-glass 
structure. It was suggested in [22] that fast quantum 
oscillations of heavy rare earth ions in a double-well 
potential lead to spin polarization of 5d-conduction band 
states, and this effect appears to be very sensitive to 
external magnetic field. Turning to the analogy with non¬ 
magnetic impurity-induced dynamic moments in HTSC 
cuprates, it should be stressed here that the localized 
moments found in present study of Lu a: Zri_ a; Bi2 are 
static and that a strong enhancement of spin polarization 
(exhibiting a saturation above 40 kOe) is induced by 
external magnetic field. 

To summarize, we have found a formation of static 
magnetic moments with p e ffps3p B per Lu-ion in 
the vicinity of non-magnetic lutetium impurities in 
the Zr-rich matrix of Lu x Zri_ a; Bi2 dodecaborides at 
low temperatures. According to our opinion, the 
strong suppression of superconductivity in Lu 2; Zri_ x Bi2 
compounds should be attributed to pair-breaking arising 
in the vicinity of these nanosize magnetic domains. 
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SUPPLEMENTARY INFORMATION 


1) 



Fig. SI: X-ray Laue pattern from the lateral cross section of a raw single crystal Luo.o74Zro.926Bi2 
grown with [100] oriented seed. The growth direction deviation from [100] is about 3 degree. 


2) three-level Schottky relation: 


C m =N 


R 2A Efe~P AEl + 2AE%e~P AE2 + (A E 1 - A£ 2 ) 2 e _/3(AiSl+AEa) 
02 (2 + e~P AEl + e~P AE2 ) 2 


P = 1/ksT 


(SI) 


3 ) 

Tab.SI: The parameters of heat capacity and magnetization analysis performed for LuzZri-^Bi^ on figs. 2 and 3. 



C 

M 

X 

gi 

g2 

Nnsm [%) 

M-eff (|+b) 

|Asat (|Ab) 

Xo (iVkOe) 

0.007 

3.1±0.1 

7.9+0.7 

0.7+0.04 

3.3 

2.88 

2.32-10 5 

0.015 

3.310.3 

7.810.2 

2.1+0.1 

3.8 

4.3 

-4.59-10" 5 

0.016 

3.010.1 

7.6+0.7 

1.8+0.2 




0.065 

3.5+0.2 

6.8+0.1 

7.210.3 




0.074 

3.2+0. 1 

7.5+0.3 

6.8+0.6 

3.3 

3.0 

3.27-10 5 



















